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Non-local electrodynamic models are developed for describing metallic surfaces 

- - - for a diffusive metal. The electric field noise at a distance zn from the surface is 
00 ■ 

' evaluated and compared with data from ion chips that show anomalous heating with 

■ a noise power decaying as Zq^. We find that high surface diffusion can account for 

j_j ■ the latter result. 

^3 : 
^ ■ 

; I. INTRODUCTION 

' A growing number of groups are using atom-chips or ion-chips to confine cold atoms, 

in anticipation of future apphcations and fundamental discoveries Progress in cold 

Ch ■ atom physics has led, however, to extreme sensitivity to electromagnetic noise from metallic 

^ ■ surfaces since the chips contain inevitably metallic components. Atom chip experiments with 

^ O^ . conventional metal structures have demonstrated the existence of magnetic fluctuations way 

above the level of blackbody radiation [2I, [sl, as predicted by a theoretical study 

^ ! Experiments with ion chips have shown "anomalous heating" ^i.e. the ions absorb energy 

CN ! faster than expected from Nyquist noise of a good metal jl, H 0, S]- The dependence 

^ i on distance zq to the surface (~ Zq^) also does not conform the expected form of RefQ, 

^ — J I proportional to z^'^. A tentative explanation are charge fluctuations in the surface of the 

^ ' metallic electrodes that enhance the electric field noise [511. One of the motivations of this 
O ' 

I work is to develop a model for this process, based on the diffusive motion of charge carriers. 
^ ' Due to diffusion, the metallic medium responds in a nonlocal way to an applied field, which 
has been well studied in the past|9|,E3- 

The topic has recently attracted interest again to 
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^ ' understand the Lifshitz-Casimir force between doped semiconductors 11 

■ In the present work we develop a method to treat non-local electrodynamics (NLED) 

corresponding to a momentum-dependent dielectric function. In particular we allow for a 
coupling between longitudinal and transverse field fluctuations via the electrons' diffusion 
constant. Our approach is valid on length scales large compared to the electronic mean free 
path. Other small scales that are not resolved here are the Thomas- Fermi (or Debye-Hiickel) 
screening length and the Fermi wavelength. Relevant scales are the distance of observation 
from a planar surface, the wavelength in the medium (skin depth) at a given frequency 
and the diffusive path length over one field period. We shall make the approximation that 
the vacuum wavelength is very large which is justified for the distances (1 — lOO/xm) and 
frequencies (100 kHz to 100 MHz) relevant for atom and ion chips. 

We find that if the charge fluctuations are confined to the surface and have a high diffusion 
constant, then the anomalous heating observed in ion micro-traps can be accounted for. The 
reasoning for such a surface layer is either (i) the short charge screening length of the metal 
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which is of the order of one atomic unit, or (ii) a distinct electron band existing at the surface. 
The latter scenario is known for a number of metals and has been recently confirmed via an 
observation of acoustic surface plasmonsfisl. 14 . 



The outline of the paper: in section II we review the method for evaluating electro- 
magnetic noise and its relation to the ion heating data. In section III we present general 
properties of the NLED method for a normal metal allowing for dissipation and diffusion. 
In section IV we apply the method to electric noise above a diffusive surface layer, while in 
section V we evaluate the more conventional case of a continuous charge. The appendices 
complete the paper with a review of a local medium scattering problem and and a nonlocal 
calculation of magnetic noise. 



II. HEATING RATE AND NOISE 



Consider a normal metal that occupies the half space at z < 0. The aim is to evaluate 
the fluctuations of the electric field E(r,t) outside the metal, due to thermal or quantum 
fluctuations in the metal. The metal is in thermal equilibrium. An efficient procedure 
Esl . EgI uses the following steps: (i) Introduce a source dipole ae at ro outside 



the metal and evaluate the wave emitted by the dipole, E*(r,t), which is the incident wave 
on the surface, (ii) Solve the scattering problem at the surface and find the refiected wave 
W{r,t). This identifies the response function aij 



-iu)t 



(1) 



(iii) The final step is to use the Fluctuation Dissipation theorem (FDT) with the interaction 
Vint = — E(r, t) ■ ae"""^*. The relevant type of FDT is related to the Golden rule that gives 
the transition rate —> 1 for a dipole 
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where we assume metal and field in thermal equilibrium at the inverse temperature l/ksT = 
13 /h (partition function Z) and sum over the final states |/). These states are actually 
polaritonic states including the excitations of the metallic medium (see, e.g., RefJlTl). 

Consider now the emitted wave from a unit charge oscillating at position r(t) = ro+ae"""^*. 
In the limit a — > the charge and current densities are 



P(r,t) 
J(r,t) 



-a ■ V5^(r - ro)e^ 



-ILjt 



-zcja5=^(r-ro)e-*^*. (3) 
The emitted wave can be found via the vector potential A{r,t) = A(r)e~*'^*, which in 
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the Lorentz gauge is 



A(r) = I J(r')^ ^rfV = ^ 



rol 



/^2^ gik-(r-ro)-i'o|2-zo 
7 — (4) 
(27r)2 t;o ^ ^ 

where Fq = (xq, yo? 2^0)5 k is a 2-dimensional vector in {x,y) directions and Vq = 
\/k^ — {uj/cY; for an outgoing wave we choose Imv^ < and Ylev^ > 0. The last hne 
of (jl]) is known as the Weyl angular spectrum, and can be proven by showing that it solves 
(V^ + ^)A(r) = — ^J(r) with outgoing boundary condition. 

The electric field in the vacuum layer < z < zq solves (V^ + ^)E(r) = 0, hence the 
incident wave has the form 

E^(r) = I (0E^(k)e^'^-^ (5) 
where kj = (k, —wq). Since E*(r) = x V x A*(r), we obtain from Eq.(j4l) 

-iki-ro ^2 -iki-ro 

E^(k) = -27rk, X k, X a = -27r[(k, ■ a)k, - —a] . (6) 

We are left then with the task of solving a scattering problem and identifying the response 
Q;jj(a;), which via Eq.(l2]) will yield the heating rate. 



III. DIFFUSIVE METAL - GENERAL PROPERTIES 

In this section we develop and apply a theory of nonlocal current-field response, equivalent 
to a g dependent dielectric function. At low frequencies we use the constitutive J(E) relation 

as 

J = aE - DVp (7) 

where D is the diffusion constant and V ■ E = 47rp. Taking the gradient of ([7]) and using 
continuity V ■ J = —dp/dt = iup we obtain 

V ■ J = Airap - DV^p = tup . (8) 

In 3D this can be Fourier transformed to yield Aira + Dq^ — iuo = 0, showing an overdamped 
plasma mode, as expected; it is also consistent with e/,(g, cj) = where the q dependent 
longitudinal response is eL^q.uj) = 1 



Defining p(r) = J ■^|^e*'^' ''p(k, 2), the solution of ([8]) has the form 



p(k, z) = p(k, 0)e"i" ^1 = Y ^^^^jj^ + P , Re t;i > . (9) 

The typical scale of is determined by the Thomas- Fermi screening length oq = {D/AnaY^'^ 
which in normal metals is ~ lA. Adopting a Drude model, we have Og ~ T^fi? / {Ame^kF) 
with kp the Fermi wavenumber. 
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Eq.Q is valid for small gradients, hence for a normal metal with a short we propose 
the following more general description. We assume that the charge has a surface component 
j{x,y)6{z) whose charge density ■y{x,y) will be determined self consistently. This surface 
layer takes care of the rapidly varying charge. In addition, we allow for a continuous charge 
p(r) that extends into the bulk z < 0. The total charge is then 

Ptot{r)=p{r)e{-z)+-f{x,y)6{z) (10) 

The divergence V ■ E involves a jump across the surface as well as a bulk term, hence 

V ■ E = (Ef * - )5(z) + (V ■ E)^(-z) = 47rp{r)e{-z) + 47r-f{x, y)5{z) (11) 

and therefore 

ET'-Ef = 47r7(x,|/) 

V ■ E = 47rp(r) z<0 (12) 

We generalize Eq.(I71) to allow for surface diffusion 

J = (aE - DVp)e{-z) - V||7(x, y)6{z) (13) 

where V|| is a gradient with components parallel to the surface. We neglect here a conductive 
surface current as the surface layer is narrow. This current would be of the form (TsaQE\^6{z) 
with Co the layer thickness and a.^ of order a. It would add to the following Eg. (1141) a term 
cr^aoVii ■ E <C (tEz- For the latter estimate we consider the typical scale ~ I/zq and get a 
ratio ~ a^k ~ ao/zQ <^ 1. Eq. f|T3l) is the simplest diffusion model that allows for enhanced 
surface diffusion relative to bulk (coefficient Dg > D) and that takes care of the broken 
isotropy at the surface. 

Charge conservation yields 

V-J = {aV DV'^ p)9{-z) - Jl^'S^z)- DsVl'^{x,y)5{z) = iuj[pe{-z) +-f{x,y)6{z)] (14) 

where J*" = ct-E™ — Ddzp{0) is the bulk current that flows into the charge layer. The 0{—z) 
terms in (1141) reproduce the bulk form ([H]) while the 6{z) terms provide the following link 
between the bulk and surface charge densities: 

JT = aET - Dd,p{0) = -DsV^^ix, y) - iu^{x, y) (15) 

In the following we consider two limiting cases: (i) The charge layer model where in the 
bulk p(r) = 0; the surface charge is then determined by (fT5l) . see Eq. (!2T|) below, (ii) The 
second case is the continuous charge model where ■y{x,y) = 0. In the latter case Eq. (fT5|) 
yields the obvious boundary condition J™ = 0. 

We proceed to derive the boundary conditions that will be used in the two models. The 
wave equation for the electric fleld is (for model (i) replace D by Dg in the following). 
Maxwell's equations are 

^ iu 

V X E = —B 

VxB = ^E + ^J (16) 

c c 
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Hence 

V X V X E = — (E + — J). (17) 



We insert the current density ([7j) and eliminate ptot by V ■ E = 47rptot that apphes to all 
z in Eq. (fT2|) . giving 

V^E + ^^e(^)E = (1 + ^) V(V ■ E) . (18) 

with e{uj) = 1 +4:Tiia/u! the usual transverse dielectric function. To estimate the correction 
Du/c^, we introduce the skin depth 5 = c/ \p2/Kauj as the length scale associated with Ohmic 
damping of the field inside the bulk [1/5 = (co'/c)Im e(u;)]. Usually, Duj/c^ = 2(ao/5)^ <C 1 
so that the main effect of the diffusive currents is via the boundary conditions. 

We have already noticed the jump in E^, Eq.Q. Due to the diffusive surface current, 
there is also a jump in the magnetic field, e.g.: 

AttD 

5» - 5-* = -^d^j (19) 

By applying the x-derivative and combining with the corresponding equation for B^, we get 
ET' - eET + ^VS7(x, y) - ^9,p(0) = . (20) 



which can also be derived by combining Eqs. fll2|15p . Finally, with the usual argument that 
■ di = J V X E (is for a contour approaching the boundary, we deduce = E^^^y. 
Other boundary conditions depend on the model for the charge layer. 



IV. CHARGE LAYER MODEL 

The charge below the metal surface is represented here as a surface charge. The relevant 
material parameters are the surface diffusion constant Dg and the bulk conductivity cr. The 
surface charge density 'y{x,y) is determined self-consistently from the boundary conditions. 
Setting the volume charge density p(r) = in Eq. lfTSl) . we find in k-space 

Recall that the wave vector k is two-dimensional, with components parallel to the metal 
surface. Eq. fl20|) yields now 

pout _ -pin ~ _ e + iDsk'^/u 

The effect of diffusion is to introduce the correction Dsk'^/uj, which for the important 
scale k = 1/5 defines the dimensionless parameter 

-t^o = — ?^ — 5 — — 71 FT • 
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Taking typical values of room temperature metal conductivities and diffusion constants in 
the bulk, and assuming ^ D, we find Dq ^ 1. We note, however, that the surface 
diffusion Dg is not well known. Furthermore, at lower temperatures Dq is significantly 
enhanced. 

To complete the boundary conditions, we recall that V ■ E = ik ■ E|| + dzEz where the 
tangential vector E|| is continuous across the interface. Considering that V ■ E = at 2; 7^ 0, 
we thus find 

a.^r* - dzEf" = 0. (24) 
Integrating the x component of (fT8|) over the surface layer leads to 



dzE:-' - dzET = (1 + '-^YKl = (1 + '-^)^^—^^k^ET , (25) 

which can also be found from the jump in By due to the diffusive surface current, Eq.f|T9 
We proceed now to solve the scattering problem. At 7^ 0, E solves 



(V^ + -)E(r) = 



C2 



2 



(V + ^e(^))E(r) = (26) 
The general form of the incident, reflected and transmitted waves is then 

E(r) = /"-^[E*(k)e^''-"+''°^ + E''(k)e^''-'-''°^] z>0 

E(r) = j ^E*(k)e^''-+- z < (27) 



where u = -y/ /c^ — ^e(a;) and the sign is chosen so that transmission decays, i.e. Re(f ) > 0. 
Boundary conditions for the z component yield 

iEi = EI + e: 

vEl = vo{El - El) (28) 

so that 

£;:(k) = ?^E:(k). (29) 

We observe that this has the same form as the Fresnel reflection coefficient in p-polarization 
(transverse magnetic), but with a nonlocal, effective permittivity e = e{k,uj) in place of the 
local value. The imaginary part of this nonlocal reflection coefficient is plotted in Figil] (thick 
solid line). This quantity determines the spectral noise power of electric fields polarized 
normal to the surface (see Section [Tll)- A clear enhancement at intermediate fc- vectors is 
seen. 

Working out the angular integration in Eq. (l27I) . the reflected field, evaluated at the source 
position, becomes 

Eliro) = a, / dke'^^-^- , (30) 
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FIG. 1: Imaginary part of TM-polarized reflection coefficient vs. wave vector k for a charge layer 
with lateral diffusion. Thick line: charge layer with lateral diffusion constant D^j the reflection 
coefficient is defined as the ratio E^i^i) / E\{}i) in Eq. ([2U|) . Thin black line: continuous charge model, 
Eq. ()49p . with bulk diffusion constant D = Dg- Thin red line (mainly constant): local calculation, 
Eq. ()A6p . k is scaled to the inverse skin depth 1/5. The imaginary part of the reflection coefficient 
is multiplied by Ime because in the local limit, Imr^ — > 2/Ime at large k. The parameter choice 
Lo5/c = 0.02 is arbitrary. 

which is the generalization of previous results obtained in the local approximation (see 
Appendix |Al). The surface response function azz{'^) can be read off here, according to 
Eq.dH). 

We now compute the behavior of Ima^zl^) short distances Zq ^ \ which determines 
the spectral strength of electric near fields. Note that the typical A;- vectors in the integral fl30|) 
are of the order of /c ~ l/^o where vq ^ k because zq is much smaller than the wavelength. 
In the regime 5 -C ^ A we have v ^ [1 — i)/6. We also neglect the Dg term in 
e - iDsk'^/uj = 1 + ^(1 - ^) since Dsk'^/Ana ~ al/ zl < 1. Hence 



A-KU 

evo — V 2 l+iDgk'^/uj 1 



evQ + V k Ania/uj 6 



(31) 



where the last terms are small, ^ ~ ^^2% ^ 1, ~ ^ ^ ^ corresponding to <^ 
zqS <^ X^. The imaginary part has then two terms 

where we used the dimensionless Dq defined in Eg. ( 1231) . Note the distinct power when 
the diffusion term dominates. 

In the regime zq S we use v ~ ^(l — i/6'^k'^) so that 

evQ — V 2 2i 



evQ + V i 6'^k'^e 

Similar calculations give 



(33) 
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FIG. 2: Electric field noise for the surface charge model, vs. observation distance zq (scaled to 
the skin depth). We plot the quantity lmazz{zo,uj) at fixed frequency and arbitrarily scaled. The 
surface diffusion constant is proportional to the dimensionless parameter Dq. Dotted line: local 
calculation by numerical integration of Eq. ijASP : solid lines and symbols: nonlocal calculations. 
For Do = 10, the symbols are found by numerically integrating the imaginary part of Eq. (l30p . the 
solid lines from the asymptotic formulas (j32p (for z > 6) and (j34p (for z < 5). Dq = 100 (thick 
solid line): numerical integration. For this plot, we arbitrarily choose uj6/c = 2tt6/X = 10~^. 



so that surface diffusion enhances the 1/zq term. These formulas are illustrated in Figj2] and 
compared to a numerical calculation based on the integral fl30|) . 

Consider next the scattering equations for the tangential field. Without loss of generality, 
we work with and have 



K + K 

Vo{El - E:) - vEl 



Ei 



1 



iK\.Et = iK,, 



1 



2vq 



l + iDgk'^/uJ evQ + V 



El 



where DgUj/c^ <^ 1 is neglected. Inserting E^. and El from 



2tt 



cPk e" 



(27r) 



^^0 



^ + 



Using e ^ e/(l + iDgk^ jio) we have for the correction to the local result ( 1A15P 



dk k- 



(35) 



(36) 



(37) 



iQ vq {v + vq) {evo + v) [evo + v{l + iDsk'^/uj)] 

For 6 zq <^ X we use v ~ ^ ^ k and ek ^ v^vD^k'^/uj. For the imaginary part we 



need the next order term which comes from v + k 
all other ones, 



1 + YZi)'i this correction is larger than 
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so that 

For zq <^ 5 we use v ~ k{l — i/5'^k'^) as the leading correction so that 

Ima = r dke~''^^- ~'^'^''^ (l-^) = -^ (40) 

To summarize all the results, we have at large distance 5 <^ zq <^ \ 

OJ I 3 ^ 52 



while at short distance zq 5 we have 



1 



Ima^^ = 3(1 + Do 

Ima,, = -^^3(1 + Do). (42) 
IQ-na Zq 



The diffusion factor Dq is seen to increase the fluctuations, in particular it introduces novel 
power laws that apply in intermediate distance ranges. 

This result may account for the anomalous heating observed in cold ion experiments 
. In fact, the surface diffusion coefficient determines a length scale ~ {Dg/ujY''^ 
for the surface charge. Our charge layer model predicts an enhanced electric field noise for 
distances z^ below 5, with a I/^q dependence, an intermediate range 

with a 1/zq dependence, while at long distances Id ^ <^ \ the noise is insensitive to the 
charge layer. Both the enhanced noise and the power law 1/zq quahtatively agree with the 
observed ion heating, provided one takes a large value Dq 3> 1 for the normalized surface 
diffusion. 

We note that regarding one point, these estimates are at variance with the previously 
invoked picture of pateh eharge fluetuafousB: the typical patch size has been thought to 
be much smaller than zq to account for the l/^g law- Here, the increase in noise sets in if 
the diffusive length -^d is larger than zq. It would be interesting to investigate further what 
length actually best represents the size of a charge patch. Another candidate could be the 
(surface) mean free path f^fp with Ds = ^"^{p/Tg [ts'. surface scattering rate]. At the low 
frequencies relevant here, this scale is indeed much smaller, f^fp = Id^/'j^Ts to- 



V. CONTINUOUS CHARGE MODEL 



In this subsection we consider the solution for the charge density Eq.Q as it is, i.e. 
without a surface charge. This is relevant if the electron density is very low so that ao is 
large compared with atomic scales, e.g. as in doped semiconductors [ill . 
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The additional equation for p(k, z) necessitates an additional boundary condition: as 
discussed below Eq. (fT5|) . the z component of the current vanishes at the surface, 

J,(k, z = 0) = cT£;f (k, 0) - Dd,p{k, 0) = . (43) 

This determines the charge density, as discussed in Eq.Q above: 

p{k,z) = ^ET{k,0)e^'^\ z<0 (44) 

and p(k, 2;) = for 2 > 0. There is a jump in the charge density, but no surface charge 
~ 6{z). 

Splitting the electric field in Eq.l fTSl) into longitudinal and transverse parts, El^t, we 
obtain 

4-7ii D 

E^{k,z) = ^-v,pMe^^^ (45) 

and 

E'^(k,^) =E^(k,0)e''^ (46) 

As expected, these fields vary on two distinct length scales. 

Since there is no surface charge, the electric field is continuous at 2; = 0: 

= * (47) 

but its derivative jumps 

a^Ef - a^Ef * = 47rp(k, 0) (48) 

as can be also seen by integrating Eq.( |T8l) across the surface. Note here that the diffusion 
term is confined to the inner part of the surface where p(k, z) is continuous; hence its 
contribution vanishes. 

The solution of the scattering problem for the z component is 

whose imaginary part is plotted in Figil] as thin solid line. It can be seen that diffusion 
in the bulk reduces the noise power on small scales (compared to the local calculation, red 
solid line), which is essentially a screening effect. 
For zq <^ A, we get a response function 



dke[l-j^^-{l~^/-y'^ (50) 



where the first two terms in the bracket correspond to a local medium (see Appendix 
For both 6 zq and 5 we find that the corrections to the local form are of order ao/zo 
and are therefore negligible. Significant changes only appear for z < Qq where the divergent 
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power laws are regularized, see, e.g. Ref.[18|. This regime is irrelevant, however, for atom 
and ion chips based on good conductors, because of the smallness of gq. 

For the E^. response, we have the usual boundary condition E!^ = E°^^. Since there is 
no surface current, the magnetic field By is continous as well, and we have dzE^^ = dzE°^*. 
Therefore 

El = l^L^JSi - '"'(^-^ '"'-"> ^"""K,^^ (51) 

vo + v [vo + v)vi evo + V + [e - l)k'^ /vi 

We find again that the corrections in \m.axx to the local theory (see Eq. flA17p ) are smaller 
than Dq by factors 0{aQ/ Zq, Zq/ X), hence they are negligible. 

We conclude that this "ideal" surface differs from the surface charge layer of the previous 
subsection. The charge layer model averages on the short scale oq and represents in some 
sense a rough surface, at which the surface charge is found in a self-consistent way. We have 
checked that the reflection coefficients found here are in agreement with the approach of 
Kliewer and Fuchs (see, e.g., Refsjo], [lo|), provided one uses the nonlocal dielectric function 
for the bulk mentioned after Eq.Q. A very similar calculation with the same results recently 
appeared in Ref . 
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VI. DISCUSSION 

Charge diffusion in the surface layer is seen to increase the fluctuations of electric flelds. 
In particular at large distances, it changes the distance dependence of the noise power from 
l/^Q into l/^Q or 1/zq, depending on the fleld polarization. This result may account for 
the anomalous heating observed in cold ion experiments In particular Refso 

and flnd power laws l/^o the electric noise with exponents between n = 3.47 ± 0.16 
and n = 3.8 ± 0.6, consistent with our prediction 1/zq from Eq. (l4T]) . The heating rate 
is 10^ — 10"^ higher than what is expected from thermal noise. This would imply enhanced 
surface diffusion (dimensionless coefficient 3> 1). These experiments also flnd a frequency 
dependence S{lj) ~ |^-o.8±o.4 ^^^^ S^u) ~ 1/u;, respectively, which is marginally consistent 
with u;"^/^ of Eq. lHTj) if the correction in Dq dominates [note the Bose factor ~ T/cj from 
Eq.(l2])]. More recent datajsj claim, however, a spectrum S{uj) ~ ^-i-4±o.4 fixed zq = 40/im 
and a noise level smaller by about two orders of magnitude compared with RefH 

The model of a charge layer represents an average over details of the charge distribution, 
and in this sense it is a model of a rough surface. It may correspond to the "patch" model 
[sl, 0] where random metallic segments contribute to the noise. The diffusion constant D 
represents correlated charge fluctuations on a scale ~ a/ D/lo which could be a patch size. 
To account for the magnitude of the noise, our parameter Dq needs to be large. A study of 
the surface diffusion is needed to identify the value of Dq. An alternative scenario leading 
to a surface charge layer is the presence of surface electronic bands in certain surfaces of 



certain metals 13 
surface plasmons 



141 ]. These surface states have led to the recent discovery of acoustic 



141 ]. In both cases, the heating rates observed in ion microtraps may be 



used as a probe of enhanced surface diffusion. 
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In a separate work we have studied surface plasmons by using similar NLED models 



191]. Furthermore, surface plasmons yield a small but finite electromagnetic noise even for 



superconductors at temperatures well below their critical temperature. 



APPENDIX A: LOCAL ELECTRODYNAMIC THEORY 

For completeness, we rederive here the results for a local €{uj), i.e. D = Dg = 0. 

Within the metal or the vacuum e{uj) is uniform so that V ■ E = 0, however at the 
boundary we have from Eq. lfTSl) (with D = 0) V[e(2;)E] = 0, where e{z) jumps from e to 
1. Integrating from z = —0 to z = +0 gives = ei?*". As usual, ^ E ■ di = yields 
E^* = E^, so that 

r-+0 

V-E = Ef*-£:f = (e-l)£:f . (Al) 
Hence there is a charge sheet at z=0 with 

V-E={e-l)ET[x,y)5{z) (A2) 

From Eq. ffTSj) . we have 

V^E + — eE = V[V ■ E] = (e - l)V[5{z)E^] . (A3) 

Integrating the z component across the boundary yields dzE°^^ —d^E^^ = 0, while integrating 
the X component yields a jump in dzE^, summarized in the following boundary conditions: 

rpin TT'Out ^ rpin rpout 

^x,y = ^x,y ' (^^z = 

dzE^^; - dzEll = (e - l)d.,,yET, dzET = dzET' ■ (A4) 

The jump in dzEr^ is actually equivalent to the continuity of the magnetic field. 

Proceeding from Eq. fl27l) we consider first the z component with boundary conditions 

eEliX) = El{k) + E:{k) 

vEiik) = vo[Elik) - Elik)] (A5) 
The solution for the reflected wave contains the Fresnel coefficient in TM-polarization, 

SO that the reflected fleld ([6]) becomes 

J {2-Kyevo + v & vo 

Terms odd in k^, ky vanish, hence the flnal form is 

E:{v,)=az dkk''-^^- . (A8) 

Jo evo + v Vo 
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As a concrete application we use the form 

e{u) = 1 + Aiiai/u (A9) 



and define the skin depth 5 = c/\/2nauj, so that for low frequencies u <^ a 



P-^^eM-^F-| (MO) 

The contribution of /c < u;/c is ~ where b = 0(1) and the radiation wavelength is 

A = 2'7ic/uj, showing an exponential decay at Zq ^ A. In the following we assume Zq ^ \ 
and consider two regimes: 



l-i 

,iiuii laiigc lias A. ~ -L / ^0 ^ i^iu / aiiu. u ~ 

(using Re{v) > 0). Expansion in 1/e yields 



(i) 6 <^ Zq <^ \: The dominant integration range has k ~ l/^o ^ euj/c and v ^ ^ 



(ii) Zq 6 'ti X: Here /c ~ I/zq ^ e^^/c so that 

= a, / dkk^- — -e-^^'^' = a,(l-') (A12) 



Hence 



u; 1 ^ 
Ima^^ = - — — ^ d <^ Zq 
Sna OZq 



^ zo<^6. (A13) 



WO" 2, 







Consider next the x component with boundary conditions ( ]A4I) 



^:+^: = El 

vo{El - El) - vEl = (6 - l)zk,Ei = (1 - h^^k,^^^El (A14) 
Eliminating i?^ yields 

ElirQ) = dkk \^v^ r - 4^ (A15) 

" Jo vo % + ^^o (i; + i;o)(ewo + f) 2c2i; + Wo^ 

This form is rather tricky for expansion, due to cancellations between the first two terms. 
We therefore rewrite these as 

v-vo 2P(e-l) v'^ + evov-vov-evi + 2k'^{e-l) evo - v 

^ 7 77 7 = 7 77 ^ = (Alb) 

V + vo [v + vq) [evo + V) [v + vq) [evo + v) evo + v 

using 2k'^{e — 1) = 2evQ — 2v'^. Hence 

ro = ia, / rf/c/c K A17 

Jo evo + v & v + Vo 
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as 



m RefH (see their Eq.(A5)). Expansion as above yields 



zo < 5 ■ (A18) 



IGvrcr Zq 



The form for 5 ^ differs from Refj^ [Eq.(A7)] by a factor 2, but it does agree with Refi 
[Eq.(20)]. 



APPENDIX B: MAGNETIC FLUCTUATIONS 



Magnetic fluctuations from a metallic surface are well studied j2|, 0] and the data is in 
good agreement with the local theory 0, 16|. It is therefore important to study the effects 
of NLED on the magnetic fluctuations and check if the electric noise enhancement when 
Ds/ D ^ 1 is still consistent with a negligible magnetic noise due to the diffusion Ds. 

Consider a magnetic moment m5^(r — ro) as a source of radiation with frequency uj. This 
source is equivalent to a current source J = V x [m5^(r — fq)], which from Eq.(j4]) yields an 
incoming electric field at < 2; < ^o, 

E,(k) = -^k, X me-^''-'-« (Bl) 



It is convenient to work with the boundary conditions of electric fields as in section IV, and 
at the end find the reflected magnetic field B^ (ro). The boundary conditions for Ez, Eq. (l22l) . 
yields 



while for the x component Eq. fl35p yields 



kyUlx) — 



-iki-ro 



V + Vq 



El - 



2Vn 



V + Vq evQ + V evQ + V 



-El 



(B3) 



where ljD/c^ <^ 1 is neglected. The reflected field is then 



Blik) = -ihEl - kyEl) = -k'- 



.V — Vq 27rm^ 



e + terms odd in k^ or ky (B4) 



00 ' " " ■ V + Vq VqC 

The odd terms vanish in the k integration, hence the magnetic response function is 



ak — 



-2vozo 



CVq Vq + V 

which is identical with the result for the local theory 0, 16 
For the reflected x-component, we have 

Blik) = -ikyEl - ivoE;) 



(B5) 



evo + V 



evo 



2vi 



- 1 



V + Vq 



<y - Vq 

'v + Vq 



-'iki-ro 



(B6) 
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up to terms odd in ky that vanish in the k integration. The change in the response function 
AB^^ = B^^ - B^^{Ds = 0) is 

Im AB,, = Im r !tl£l!fhL -""t'Vt ■ 

Jo cvq v + vq [evQ + v[l + iDsk^/uj))[evQ + v) 

For 5 <^ zq <^ \ we have Wq ~ /c, f ~ ^ ^ /c ~ l/^o, ev^ ^ hence 

ImAiS,, ^^-^ (B8) 
ca ozq 

which is smaller then the local theory result jl, [l^ ^/c^o by a factor of a^Ds/S'^D ^ 
10-^2£)^/£) ^i^h a typical 5 ^ 100/i. For 5 > we have v ~ A;(l — i/S^k"^), hence 

ImAS,, ^— (B9) 
ca d^z^ 

which is small compared with the local theory result l/c^o^^ by a factor of DquS"^ /az^ ~ 
10~^^-Do with 6 ~ zq. We conclude that even if Dq is sufficiently large to account for the 
observed electric noise, e.g. Dq ~ 10^ — 10"^, the effect on the magnetic noise is negligible 
and therefore our NLED model of a charge layer is consistent with the magnetic noise data 

i,a. 
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